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Abstract 

We study the convergence of a drift implicit scheme for one-dimensional SDEs 
that was considered by Alfonsi [T] for the Cox-Ingersoll-Ross (CIR) process. Under 
general conditions, we obtain a strong convergence of order 1. In the CIR case, 
Dereich, Neuenkirch and Szpruch [2] have shown recently a strong convergence of 
order 1/2 for this scheme. Here, we obtain a strong convergence of order 1 under 
more restrictive assumptions on the CIR parameters. 
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This paper analyses the strong convergence error of a discretization scheme for the 
Cox-Ingersoll-Ross (CIR) process and complements a recent paper by Dereich, Neuenkirch 
and Szpruch [2] . The CIR process, which is widely used in financial modelling, follows the 
SDE: 

dX t = (a- kX t )dt + a^X t dW u X = x. (1) 

Here, W denotes a standard Brownian motion, a > 0, fcGl, a > and x > 0. This SDE 
has a unique strong solution that is nonnegative. It is even positive when a 2 < 2a and 
x > 0, which we assume in this paper. It is well-known that the usual Euler-Maruyama 
scheme is not defined for ([I]). Different ad- hoc discretization schemes have thus been 
proposed in the literature (see references in [2]). Here, we focus on a drift implicit scheme 
that has been proposed in Alfonsi pp. We consider a time horizon T > and a regular 
time grid: 

kT 

tk = — , < k < n. 

n 



By Ito's formula, Y t = \fX~t satisfies : 



dY 



a - a 2 /A k 



a 



^ _ Yt \dt+-dW t , Y = y/E. 



(2) 



We consider the following drift implicit Euler scheme 

'a-a 2 /A 



Y = y/j,Y t = Y tk + 



2Y t 



\y\ (t-t k ) + ^(W t -W th ),te(t k ,t k+1 }. (3) 



The equation (J3J) is a quadratic equation that has a unique positive equation: 



Y tk + §(Wt - W tk ) + J(Y tk + |(W f -W tk )) + 2 (1 + |(t - i*)) a - ^ (t - t fe ) 



2(l + |(t-* fc )) 



provided that the time-step is small enough (T/n < 2/ max(— fc, 0) with the convention 
2/0 = +oo). Last, we set X t = {Y t ) 2 , t G ifc+i]. It is shown in p] that this scheme has 
uniformly bounded moments. We recall now the main result of Dereich, Neuenkirch and 
Szpruch [2] that gives a strong error convergence of order 1/2. 

Theorem 1. Let x > 0, 2a > a 2 and T > 0. Then, for all p E [1, t/iere zs a constant 
K p > sitc/i i/iai /or any n > ? max(— fc, 0), 



E 



max \X t — X t \ p 
te[o,T] 



i/p 



< K r 



Let us remark that, contrary to [2j, we do not consider a linear interpolation between t k 
and tk+i here for X t . This removes the logarithm term of Theorem 1.1 in [2]. 

The strong convergence rate of X is studied numerically in Alfonsi (pQ, Figure 2). This 
numerical study shows that the strong convergence rate depends on the parameters a 2 and 
a. When a 2 /a is small enough, a strong convergence of order 1 is observed. The scope of 
the paper is to prove the following result. 

Theorem 2. Let x > 0, a > a 2 and T > 0. Then, for all p e [1, -^2), there is a constant 
K p > such that for any n > ^ max(— k, 0), 



E 



max \X t 

te[o,T] 



XA p 



i/p 



n 



Thus, we get a strong convergence of order 1 under more restrictive conditions on a 2 /a. 
Both theorems are complementary and are compatible with the numerical study of PQ, 
which indicates that the strong convergence order downgrades as long as a 2 /a increases. 



The paper is structured as follows. We first prove that ( E 



maXi e[0i T] \ Y t - Y t \ p 



i/p 



< 



K p — under a general framework for Y and Y that extends (j2j) and ([3]). Then, we deduce 
Theorem [2] from this result. Also, we construct an analogous drift implicit scheme for 
general one-dimensional diffusion, and get a strong convergence of order one under suitable 
assumptions on the coefficients. This scheme has the advantage to be naturally defined in 
the diffusion domain like Ml for the CIR case. 
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A general framework for Y and Y 



Let c G [—00, +00), I = (c, +00) and d G I. We consider in this section the following SDE 
defined on I = (c, +00): 

= f(Y t )dt + >ydW u t>0,Y = yeI, (4) 
with 7 > 0. We make the following monotonicity assumption on the drift coefficient /: 

/ : I -> R is C 2 , such that 3k G R,\/y,y' el, y< y', f{y') - f{y) < K {y' - y). (5) 
Besides, we assume 



v(x) 



exp 



r 



f(£)d£ I dzdy satisfies lim v(x) 

J x— >c+ 



-00. 



(6) 



The Feller's test (see e.g. Theorem 5.29 p. 348 in [5]) ensures that Y never reaches c nor 
+00 by <^j, and the SDE (jl]) admits a unique strong solution on /. 

Let us now define the drift implicit scheme. Let us first observe that for h > such 
that nh < 1, y i-> y — hf(y) is a bijection from / to R. Indeed, it is continuous and we 
have from (J4J): 

y<y',y'-y- h(f(y') - f(y)) > (1 - - y). 

This shows the claim for c = —00. For c > —00, we first remark that lim c+ / exists 
from §5§, and is necessarily equal to +00 from (jH]). Thus, for n such that kT/u < 1, the 
following drift implicit Euler scheme is well defined 

Y Q = y, Y t = Y tk + f(Y t )(t-t k )+^W t -W tk ), te(t k ,t k+1 ], 0<k<n-l, (7) 

and satisfies Y t G I, for any t G [0,T]. From a computational point of view, let us 
remark here that in cases where Y tk+1 cannot be solved explicitly like in the CIR case, 
Y tk+1 can still be quickly computed from Y tk and W tk+1 — W tk thanks to the monotonicity 
of y 1 — y y — (T/n)f(y) by using for example a dichotomic search. 

The drift implicit Euler scheme (also known as backward Euler scheme) has been studied 
by Higham, Mao and Stuart jl] for SDEs on R d with a Lipschitz condition on the diffusion 
coefficient and a monotonicity condition on the drift coefficient that extends (jHJ). They 
show a strong convergence of order 1/2 in this general setting. 



Proposition 3. Let p > 1 and n > 2nT . Let us assume that 



E 



T 



\f'(Y u )f(Y u ) + ^f"(Y u )\du 



< 00 and E 



T \ P/2' 

(f'(Y u )fduj 



< 00. 



Then, there is a constant K p > such that: 



E 



max \Yf — YA 
te[o,r] 



i/p 



T 

<K P -. 

n 
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Before proving this result, let us recall that the same result holds for the usual (drift 
explicit) Euler-Maruyama scheme when / = R (i.e. c = — oo), under some regularity as- 
sumption on /. Said differently, the Euler-Maruyama scheme (Y tk = Y tk + f(Y tk )T/n + 
liWtk+i ~ ^tfc)) coincides with the Milstein scheme when the diffusion coefficient is con- 
stant, and its order of strong convergence is thus equal to one. The main advantage of the 
drift implicit scheme is that it is well defined when c > — oo while the Euler-Maruyama is 
not, since the Brownian increment may lead outside I. 

Proof. We may assume without loss of generality that k > 0. For t G [0,T], we set 
e t = Y t - Y t . From ©, there is A < k, such that f(Y t ) - f(Y t ) = f3 t e t . For < k < n - 1, 
we have 



et k+1 = e tk + [f(Y tk+1 )-f(Y tk+1 )]- + 



71 



■k+l 



f(Y tk+1 ) - f(Y s )ds, 



and then, by using Ito's formula: 



T 

) e tk+1 =e tk + 



'k+l 



(u-t k )[f\Y u )f{Y u ) + 1 -f\Y u )}du+ 1 



k+l 



tk 



{u-t k )f(Y u )dW u . 
(9) 



For u G [0,T], we denote by rj(u) the integer such that t^u) < u < tu u )+i- We set Il = 1, 
n fc = nti(l - ^J), h = U k e tk , fl k = n fe /(l - nT/nf and 



M t = / {l-KT/nf^iu-t^fiY^dWu. 
Jo 

Let us remark that IT^. > 0, IL; > 1 and 11^ is nondecreasing with respect to k. By 
multiplying equation (j5J) by we get 



e k+ i = e k + I1 A 



•k+l 



r 



(u - t k )[f'(Y u )f(Y u ) + -L-f"(Y u )]du + I (u - t k )-ff'(Y u )dW u ) . 



■k + l 



tk 



Then, we obtain e k = f Q k U v(u) (u - t v{u) )[f'(Y u )f(Y u ) + ^f"(Y u )]du + £L"o IT Z (M <;+1 - M k ) 



by summing over k and finally get 



e tk 



^ ^% - t n{u) )[f(Y u )f(Y u ) + (Y u )]du W +1 -M tl ). (10) 

Ok k l=Q 



Since < exp(2x) for x G [0, 1/2], we have 



0<l<k<n, 0< 



n. 



n. 



< cxi> ( 2(k-l)n— ) < <<xp{2K.7 : ). 



Tl k (1 - nT/nf- 1 n fe 

On the other hand, an Abel transformation gives Yli=o n^(M t — M t; ) = U. k _iM tk + 
Eti (fi,_i - tli)M tl and thus 



fc-i 



1=0 



k-l 



< n*_i|M t J +y)(n,-n,_i)|M tl | < 2ft* max \M tk \ 



since Ilk is nondecreasing. From (fTOl) and ^ = ( 1 _J r / n - } k < exp(2ftT), we get 

\e tk \ < exp(2*;T) (- f \f'(Y u )f(Y u ) + ^f\Y u )\du + 2 max \M tl \) . 

\n Jq Z o<i<k l 

Since the right hand side is nondecreasing with respect to k, we can replace the left hand 
side by max <Kfc |e t , |. Burkholder-Davis-Gundy inequality gives that 



E 



max \M t , 

0<l<n 



< C p j p (T/n) p E 



T \ P/2' 

(f(Y u )) 2 du) 



since < (1 - nT/n)^ < 1. Thus, there is a positive constant K depending on k, T and 
p such that: 



E 



max \e tl \ 

0<Kn 



< K 



T 



n 



E 



+7 P E 



\f(Y u )f(Y u ) + ^f'(Y u )\du 



T \ P/2' 

(f(Y u )) 2 du) 



It remains to show the analogous upper bound for E[max tg [o,T] |e*| p ]- Similarly to we 
have for t G [tfc,tfc+i]: 

(l-/3 t (t-t fc ))e t = e tfe + / («-i t )[/'(i;)/ft) + ^-m)]^ + 7/ (u - t k )f'(Y u )dW u . 



Since (1 - (3 t (t - t k )) > 1/2, we get: 



max \e t \ < 2 e t , H 

te[tfc,t fe+ i] V n 



+7 max 

<G[tfc,tfc+i] 



t'k 



\f(Y u )f(Y u ) + ^f"(Y u )\du 



{u-t k )f(Y u )dW u 



tk 



and thus 



p / pT 



max |e t | p < 2 P 3 P " 1 f max |ej p + f-^ ( 

;e[0,T] V,0<fe<n \ n J \ 

u-t n(u) )f{Y u )dW u 



\f(Y u )f(Y u ) + ^f(Y u )\du 



+7 P max 

0<s<i<T 



Since 



< 2 p 



r o (n-t v{u) )f\Y u )dW u P +\j;(u-t ri{u) )f'(Y u )dW u \ p ), 



we conclude by using once again Burkholder-Davis-Gundy inequality, fjTT|) and ([8]). □ 
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Application to the CIR process 



+ >, f(y) = s =S i± ~ \V and 7 = a/2. When 



For the CIR case, we have c = (i.e. I - ^ < f <<> - 

2a > a 2 , we can check that both and (jSJ) are satisfied. By Jensen inequality, 
if we have 

' E[|/'(y u )/(y u )r + |/"(y u )r + < oo 



holds 



:i2i 



The moments of the CIR process can be uniformly bounded on [0, T] under the following 
condition (see [2] equation (7)): 



sup E[X?] < oo for q > 
te[o,T] °" 



Condition ( 1T2"j) will hold as soon as sup tg r T i E[F t 



-(4V3p)i 



is satisfied when a 2 < a and p < t^, and we have ( E 



3(T 2 ' 



max tG[0 ,T] \Y t - Y t \ p 



•(2V|p)! 



(13) 



< oo. This 



i/p 



< K, 



T 



Let £ > such that 



From now on, we assume that a 2 < a and consider 1 < p < 
p(l + s) < Since X t - X t = (Y t - Y t )(Y t + Y t ), we have by Holder's inequality: 



E 



max \X t — XA 
te[o,T] 



< E 



max 
te[o,T] 



|y t _yjp(l+e) 



E 



max |F t + yj p " 

te[o,T] 



P(l+e) 



The moment boundedness of Y is checked in [I] and [2], and the second expectation is thus 
finite. Proposition E] gives Theorem |5J 



Application to dX t = (a - kX t )dt + crX?dW u with 1/2 < a < 1 

We consider this SDE starting from X = x > with parameters a > 0, G R and a > 0. 
This SDE is known to have a unique strong positive solution X, which can be checked 
easily by Feller's test for explosions. We set 

Y t = X]- a . 

It is defined on I = R+ and satisfies (jH) with 

f{y) = (!-«) ^T 1 ^ - % - "y^ 1 ) with 7 = a{l - a). 

Since a > and > 1, / is decreasing on (0, e), for e > small enough. It is also clearly 
Lipschitz on [e, +oo), and (j5J) is thus satisfied. Also, we check easily that ([6]) holds. The 
drift implicit scheme (Y t ,t G [0, T}) given by (JTj) is thus well defined for large n and we set: 

X t = (Y t )T^. 

To apply Proposition [31 it is enough to check that (fl2|) holds. To do so, we have the 
following lemma. 
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Lemma 4. We have: Wq G R, sup tg [ 0T ] E[Xf] < oo. 

Proof. For g > 0, it is well known that we even have E[max tG [o,T] Xf\ < oo from the 
sublinear growth of the SDE coefficients (see e.g. Karatzas and Shreve [5], p 306). Let 



q < 0. We set Z t = X 2 t {1 a) and have: 



dZ t = b(Z t )dt + 2(1 - a)ay/Z t dW t , with b(z) = 2(1 - a 



1-2" , o / 1 
a/2 2(l- Q ) _ ^ _|_ (J q; 



Since lim 2 _ >0 + = +oo and b is Lipschitz on [e, +oo) for any e > 0, we can find for 
any M > a constant /ca/ € R such that b(z) > M — kuz for all z > 0. We consider then 
the following CIR process: 



der = (m - fc M er + 2(1 - «w& M ^ $ 



2(l-a) 



From a comparison theorem (Proposition 2.18, p 293 in [5]) we get that Vt > 0, Z t > £. 
and thus sup te r 0T i E[Z t g ] < sup ie r 0T i E[(£ 4 M ) 9 ]. We conclude by using ( TT3"j) and taking M is 



arbitrary large. 



□ 



We can then apply Proposition [3] and get, for any p > 1 and n large enough, the exis- 

i/p 



tence of a constant K p > such that ( E 



max t£ [ ,T] \Y t - Y t \ 



we get E[max 4e [ ,T] Yf] < °°- We have X = (Y^) 1 -" and 



y 1 - 



1 — a 



z 1 ~ a dz 



< 



a 



< K p ^. In particular, 



\y-yKyvyY 



The Cauchy-Schwarz inequality leads then to 



E 



max \X t — X t \ p 
te[o,T] 



< 



-E 



a 



max \Y t -Y t \ 2p 
te[o,T] 



■2p 



E 



max(F t VY t )* 



2 pa. 



~ T 
<K P -. 

n 



Strong convergence towards X in a general framework 

Let us now consider a one-dimensional SDE with Lipschitz coefficients b, a : R — > R: 

dX t = 6(X t )dt + <7(X t )dW*, X = x. 

We will consider the Lamperti transformation of this SDE. We assume that there exist 
< a < a such that a < a(x) < a, so that 



cr(z 



-dz is bijective on 



Lipschitz and such that ip 1 is Lipschitz. Besides, we assume that a G C 1 and that 
/ = — y) o ip- 1 satisfies ©, © and: 



3K>0,q> 0,Vy G R, + \f"(y)\ < C{1 + |y|*). 
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Then Y t = ip(X t ) satisfies dY t = f(Y t )dt + dW t . The Lipschitz assumption on the coeffi- 
cients b and a ensures the boundedness of moments of X and thus of Y. The condition (JSJ) 
is thus satisfied and the conclusion of Proposition [3] holds. Then, defining Y by (E]) and 
setting X t = y- l (Y t ) for t e [0,T], we get that: 



3K P > 0, E 



max \X t -X t \ p 
te[o,T] 



l/p rp 

< K, 



p 
n 



Let us mention that the same result holds under suitable conditions on / for the scheme 
X t = ip~ l (Y t ), where Y denotes the Euler-Maruyama scheme dY t = f(Y t ,)dt + dW t . 
The weak convergence of this scheme has been studied by Detemple, Garcia and Rindis- 
bacher [3]. 

Remark 5. Let 7 > 0, <f-y(x) = ^^(x) and f 7 (y) = 7/(27/7)- Then, Y"/ = <£> 7 (Jf t ) solves 
dYl = f^(Y t ')dt + jdWt- The associated drift implicit scheme 

f ' = ^ 7 (X ), % = % h + f y (Y;)(t-t h ) + 'y(W t -W th ), te(t k ,t k+1 ], 0<k<n-l, 

clearly satisfies Y t ' = jYt. Thus, X t = ip^iY^): the scheme X is unchanged when the 
transformation between X and Y is multiplied by a positive constant. 
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project and of the "Chaire Risques Financiers" of Fondation du Risque. 
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